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[t] In several biomes, including croplands, wooded savannas, and tropical forests, many 
small fires occur each year that are well below the detection limit of the current generation 
of global burned area products derived from moderate resolution surface reflectance 
imagery. Although these fires often generate thermal anomalies that can be detected by 
satellites, their contributions to burned area and carbon fluxes have not been systematically 
quantified across different regions and continents. Here we developed a preliminary 
method for combining 1-km thermal anomalies (active fires) and 500 m burned area 
observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) to 
estimate the influence of these fires. In our approach, we calculated the number of 
active fires inside and outside of 500 m bum scars derived from reflectance data. 

We estimated small fire burned area by computing the difference normalized bum ratio 
(dNBR) for these two sets of active fires and then combining these observations with 
other information. In a final step, we used the Global Fire Emissions Database version 3 
(GFED3) biogeochemical model to estimate the impact of these fires on biomass burning 
emissions. We found that the spatial distribution of active fires and 500 m burned areas 
were in close agreement in ecosystems that experience large fires, including savannas 
across southern Africa and Australia and boreal forests in North America and Eurasia. 

In other areas, however, we observed many active fires outside of burned area perimeters. 

Fire radiative power was lower for this class of active fires. Small fires substantially 
increased burned area in several continental-scale regions, including Equatorial Asia 
(157%), Central America (143%), and Southeast Asia (90%) during 2001-2010. Globally, 
accounting for small fires increased total burned area by approximately by 35%, from 
345 Mha/yr to 464 Mha/yr. A formal quantification of uncertainties was not possible, 
but sensitivity analyses of key model parameters caused estimates of global burned area 
increases from small fires to vary between 24% and 54%. Biomass burning carbon 
emissions increased by 35% at a global scale when small fires were included in GFED3, 
from 1.9 Pg C/yr to 2.5 Pg C/yr. The contribution of tropical forest fires to year-to-year 
variability in carbon fluxes increased because small fires amplified emissions from Central 
America, South America and Southeast Asia — regions where drought stress and burned 
area varied considerably from year to year in response to El Nino-Southem Oscillation 
and other climate modes. 
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1. Introduction 

[ 2 ] Quantitative estimates of contemporary burned area are 
needed to address many ecological and Earth science chal- 
lenges, including understanding disturbance-driven changes 
in ecosystem function and carbon fluxes [Luo and Weng, 
2011], the impacts of land use on atmospheric chemistry 
and human health [. Johnston et al., 2012], and the strength of 
climate-carbon cycle feedbacks during the twenty-first cen- 
tury [e.g., Golding and Betts, 2008; Kloster et al., 2012], 
Global-scale burned area products are generated using 
change detection algorithms applied to time series of surface 
reflectance from moderate resolution spectrometers, typi- 
cally at 500 m to 1 km resolution, and require calibration and 
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validation with higher resolution imagery [ Roy et ah, 1999; 
Stroppiana et ah, 2003; Silva et ah, 2003; Tcinsey et al., 
2004; Simon et ah, 2004; Roy et al., 2005; Loboda et al., 
2007; Tansey et al., 2008; Giglio et al., 2009; Roy and 
Boschetti, 2009]. Burned area time series have been used as 
drivers of biogeochemical and fuel load models, often fol- 
lowing the conceptual framework developed by Seiler and 
Crutzen [1980], to estimate carbon emissions and other 
trace gas and aerosol fluxes [Scholes et ah, 1996; Barbosa 
et al., 1999; Amiro et al., 2001; Hoelzemann et al., 2004; 
Korontzi et al. , 2004; Mouillot and Field, 2005; Wiedinmyer 
et ah, 2006; Ito et al., 2008; Schultz et al., 2008; Reid et al., 
2009; van der Werfetal., 2010; French etal., 2011; Turetsky 
et ah, 2011], With this approach, information on fuels and 
combustion completeness must come from other sources, 
including field observations, satellite maps of vegetation 
cover, or model estimates. 

[ 3 ] Studies evaluating contemporary fire trends often make 
use of another complementary source of information: mea- 
surements of active fires derived from visible and infrared 
radiation emitted directly by fire fronts during the time of a 
satellite overpass [Dozier, 1981; Matson and Dozier, 1981]. 
Several commonly used active fire products rely on anoma- 
lously high brightness temperatures observed in a middle 
infrared spectral band (near 4 pm) to detect fires, and use 
additional information from visible, near infrared, and ther- 
mal infrared bands to separate fires from sun glint, clouds, 
and high land surface temperatures [Giglio et ah, 2000; 
Giglio et ah, 2003; Xu et ah, 2010]. Active fires by them- 
selves allow for the detection of the presence or absence of 
fire in a given pixel during an individual satellite overpass or 
sampling interval. When aggregated temporally and spa- 
tially, active fires allow for the quantification of temporal and 
spatial patterns [Prins and Menzel, 1992; Kasischke and 
French, 1995; Scholes et ah, 1996; Giglio et ah, 2006a; 
Giglio, 2007; Chuvieco etal., 2008; Zhang and Kondragunta, 
2008; Roberts et ah, 2009], multiyear trends [Prins and 
Menzel, 1994; Aragao and Shimabukuro, 2010; Arino et ah, 
2012], relationships between fires and climate drivers 
[Spessa et ah, 2005; van der Werf et ah, 2008], and investi- 
gations of socioeconomic and cultural factors that influence 
fire use in different countries [Mollicone et ah, 2006; 
Bowman et ah, 2011; Lin et ah, 2012], 

[ 4 ] Additional information on fire radiative power (FRP), 
the rate of electromagnetic energy released by fire, is avail- 
able for some satellite spectroradiometers whose detectors 
allow for accurate quantification of fire temperatures. Since 
the energy content of dry biomass does not vary consider- 
ably across different plant species, ecosystem type, or fuel 
type [Pooter and Villar, 1997; Chapin et ah, 2002], FRP 
measurements may enable direct estimates of fuel con- 
sumption, if they can be integrated through time to obtain 
estimates of fire radiative energy (FRE) [Ellicott et ah, 2009; 
Vermote et ah, 2009; Freeborn et ah, 2011], This relatively 
direct approach for estimating emissions may shortcut some of 
the uncertainties associated with biogeochemical model esti- 
mates of fuel loads and combustion completeness [Kaufman 
et ah, 1998; Wooster et ah, 2003], yet presents its own chal- 
lenges [Kaiser et ah, 2012], 

[5] Burned area and active fire/FRP products both have 
strengths and weaknesses with respect to their use in devel- 
oping a global time series of fire carbon emissions [Roberts 


et ah, 2011]. Burned area is likely to be less sensitive to 
variability in cloud cover and smoke in many regions since 
change detection algorithms can draw upon multiple satellite 
overpasses to obtain cloud-free scenes before and after a fire 
event. Burned area also provides direct estimate of final fire 
size and thus the total area experiencing emissions losses, 
enabling immediate use with models following the Seiler and 
Crutzen [1980] paradigm. Contemporary global burned area 
products have a limitation, however, because they are made 
using moderate resolution surface reflectance imagery, with 
pixel areas varying between 21 ha and approximately 100 ha 
[Roy et ah, 2005; Loboda et ah, 2007; Tansey et ah, 2008; 
Giglio et ah, 2009]. In some biomes there is accumulating 
evidence that many individual bum scars have sizes that are 
below the detection limit of these products [Eva and Lambin, 
1998; Laris, 2005; McCarty et ah, 2009; Rov and Boschetti, 
2009], 

[6] Small fires often occur in agricultural settings where the 
size of an individual field limits the fire spread [McCarty et ah , 
2009], where the use of prescribed fires is important for eco- 
system management [e.g., Covington et ah, 1997], where 
wildland fires are suppressed near human settlements 
[DeWilde and Chapin, 2006], in ecosystems with heteroge- 
neous patches of land cover that limit the continuity of dry 
fuels [Eva and Lambin, 1998; Laris, 2005], and in places and 
times where fuel moisture and atmospheric conditions do not 
sustain high fire spread rates [Van Wagner, 1977; Hely et ah, 
2001]. In the U.S., for example, agricultural field sizes vary 
by region, with an average of approximately 16 ha in many 
southern and eastern regions [McCarty et ah , 2009] . Assuming 
that agricultural fires bum individual fields, this places a 
regionally varying upper bound on fire sizes detectable from 
satellite imagery. Developing calibration or validation esti- 
mates of burned area associated with many of the agricultural 
and prescribed fires described above for remote sensing- 
derived products is made difficult because few state and 
national programs systematically quantify this type of burned 
area with high levels of accuracy. Further, since several burned 
area detection algorithms require that a substantial fraction of 
an individual pixel’s area undergo burning for successful 
attribution (to avoid commission errors from other forms of 
land cover change), detection of these sub-pixel fires can be 
difficult [Roy and Landmann, 2005; Roy and Boschetti, 2009]. 
As a result, the contribution of small fires to continental and 
global-scale scale patterns of burned area and carbon emis- 
sions remains poorly understood. 

[7] Active fire products have unique properties that may 
provide new information about small fire contributions to 
burned area and emissions. Specifically, these products have 
the advantage of enabling the detection of fires that are con- 
siderably smaller than the spatial resolution of an individual 
moderate resolution surface reflectance pixel because of the 
strong nonlinearity in radiative power as a function of fire 
temperature [Giglio et ah, 2003], For example, the 1-km 
MOD1S active fire product has a 50% detection efficacy for 
fires in different vegetation types that span approximately 
25 to 70 actively burning 30 m Advanced Spacebome Ther- 
mal Emission and Reflection Radiometer (ASTER) pixels 
[Morisette et ah, 2005; Csiszar et ah, 2006; Schroeder et ah, 
2008]. This area is equivalent to approximately 10% to 
30% of a 21.5 ha nadir MOD1S surface reflectance pixel. 
Given that actively burning fire fronts are not expected to 
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completely cover the area of many ASTER fire pixels, the 
detection limit for MOD1S probably is even higher — and 
more consistent with the area estimates from Giglio et al. 
[2003] based on spectroradiometer sensitivities. 

[8] The relatively high sensitivity of the MOD1S active 
fire product to small flaming fronts within each pixel makes 
it potentially useful for quantifying the spatial and temporal 
distribution of small fires at a global scale. The challenge, 
however, from the perspective of developing a biomass 
burning emissions time series is that this product provides 
only a relative measure of fire activity. Integration of FRP to 
FRE from MODIS, for example, requires accounting for 
gaps created by clouds and orbital geometry, estimating the 
diurnal cycle of FRP during periods in between satellite 
overpasses, adjusting for variable detection efficiencies as a 
function of scan angle, and removing the effects from com- 
mission errors [Freeborn et al., 2009; Freeborn et al, 2011]. 
An additional challenge with respect to the use of FRP 
involves uncertainties related to the dissipation of fire energy 
in sensible, latent, and radiant forms. Radiant energy is only a 
small fraction of the immediate energy release from a fire 
[ Wooster et al. , 2005] , and important uncertainties remain with 
respect to understanding how this fraction varies between 
different types of fires (i.e., flaming versus smoldering) and 
levels of fuel moisture [. Freeborn et al, 201 1], 

[ 9 ] Here we combined information from active fires and a 
500 m burned area product to estimate the contribution of 
small fires to global biomass burning carbon emissions. In our 
approach, we tracked the difference normalized bum ratio 
(dNBR) [Lopez Garcia and Caselles, 1991; van Wagtendonk 
et al, 2004] for all active fire pixels inside and outside of 
burned areas in the 2001-2010 MODIS record, separately 
compiling statistics as a function of vegetation type, conti- 
nental-scale region, seasonal period, and year. These data were 
then used to estimate the burned area associated with active 
fires located outside of the Giglio et al [2009] burned area 
product. We hereafter refer to this set of burned area as the 
‘small fire’ burned area because sub-500 m fires are likely to 
be a primary cause of active fires residing outside of the 
burned area product in many biomes, including agriculture, 
wooded savannas, and tropical forests. We evaluated our 
approach using FRP observations and several higher spatial 
resolution burned area products from North America, Europe, 
and Africa. To quantify uncertainties, we performed five sen- 
sitivity studies in which we modified key aspects of our 
approach or used different remote sensing data products to 
estimate small fire burned area. In a final step we investigate 
the impacts of including the burned area from these small fires 
on global fire emissions using the biogeochemical model from 
the Global Fire Emissions Database version 3 (GFED3) [van 
der Werf et al, 2010], We show that accounting for small 
fires considerably increased the total amount and interannual 
variability of carbon emissions from tropical forest regions. 

2. Methods 
2.1. Data Products 

[ 10 ] We used the MCD64A1 burned area product [Giglio 
et al, 2009] and the Collection 5 Level 3 daily composition 
thermal anomaly/fire product MOD14A1 (from NASA’s 
Terra satellite) [Giglio et al, 2003; Morisette et al, 2005; 
Csiszar et al, 2006; Giglio et al, 2006a; Schroeder et al. 


2008; Giglio, 2010] in our primary analysis. The Aqua daily 
composition thermal anomaly/fire product (MYD14A1) was 
used in sensitivity tests described below. The MCD64A1 
burned area product has a 500 m spatial resolution and 
monthly temporal resolution, although each file specifies the 
individual day of burning for every burned area pixel and a 
temporal uncertainty range for its bum date. The spatial and 
temporal coverage of the MCD64A1 burned area reported 
here is significantly greater than that used by Giglio et al. 
[2010] or subsequently by van der Werf et al. [2010] in 
GFED3: we used an updated version that included more than 
98% of all MODIS land tiles during 2001-2010. Compared 
with GFED3, this updated set of burned area had the same 
global mean (345 Mha/yr during 2001-2010) but different 
regional patterns, with decreases in Central, Southeast, and 
Equatorial Asia, and increases in Central America and Northern 
Hemisphere South America. The two thermal anomaly (active 
fire) products from Terra and Aqua have a 1-km spatial reso- 
lution, report the presence or absence of fire each day during 
one or more daytime and nighttime overpasses, and include 
additional information about the confidence level of the 
retrieval, cloud or water cover, scan angle, and fire radiative 
power. Both the burned area and active fire products are defined 
on the MODIS sinusoidal grid. 

[ 11 ] We excluded pixels with persistent non-fire hot spots 
from the MODIS daily composition thermal anomaly/fire 
products (MOD14A1 and MYD14A1) using a location data- 
base developed by Csiszar et al. [2006], gas flare pixels 
identified using the NOAA Global Gas Flare Estimates data- 
base [Elvidge et al, 2009] (Google kmz files downloaded 
from www.ngdc.noaa.gov/dmsp/interest/gas_flares.html on 
9 Mar 20 1 0), and volcano locations from the NOAA National 
Geophysical Data Center (NGDC) Volcano Data and Infor- 
mation (Significant Volcanic Eruptions Database down- 
loaded from www.ngdc.noaa.gov/hazard/volcano.shtml on 
23 September 2011). 

[ 12 ] As described below, we calculated dNBR for each 
active fire detection in the Terra time series using surface 
reflectance observations from the Vegetation Indices 16-day L3 
(MOD13A1) collection 5 product [Huete et al, 2002]. We 
chose Terra over Aqua for the surface reflectance time series 
because of its longer temporal coverage and we chose the 
1 6-day vegetation indices (MOD 1 3 A 1 ) product over the 8-day 
surface reflectance (MOD09A1) product because of the 
additional cloud screening provided by maximum compos- 
iting in the vegetation indices product. NBR is defined as 
(NIR-SWIR)/(N1R + SWIR) following van Wagtendonk 
et al. [2004] where NIR represents surface reflectance from 
a near infrared band and SWIR represents surface reflectance 
from a shortwave infrared band. For NIR we used band 2 and 
for SWIR we used band 7 from MOD13A1. Terra MODIS 
sensor degradation that impacts collection 5 surface reflec- 
tance estimates in visible wavelengths [Wang et al, 2012] 
does not impact the longer wavelength bands used in this 
study. We computed the difference in NBR (dNBR) using the 
two NBR values immediately before and immediately after 
each active fire detection (and thus excluding the 16-day 
interval containing the day of the active fire detection to limit 
complications from the evolving burned area and smoke aero- 
sols). The MOD13A1 product has a 500 m spatial resolution, 
and so NBR values for each set of 4 pixels corresponding to 
a single 1-km active fire pixel were averaged together for the 
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pre-fire and post-fire periods prior to computing dNBR. We 
used only surface reflectance pixels that had the highest quality 
retrievals — this included the set where the pixel reliability flag 
was equal to 0 (use with confidence) and the quality flag was 
equal to 00 (good quality). As described below, we also con- 
ducted a sensitivity test in which we computed dNBR for 2010 
using the Surface Reflectance Bands 1-7 (MOD09A1) collec- 
tion 5 product. 

2.2. Small Fire Estimation 

[ 13 ] To estimate monthly burned area in each 0.25° x 
0.25° land grid cell (approximately 28 km x 28 km at the 
equator), we used the following equations: 

BA tota ,(i, t, v) = BA MCD6 4 Ai{i,t,v) + BA sf (i,t,v), (1) 


such as the rate of movement of the fire front [e.g., Giglio 
et al., 2006b], a in equation (2) is not directly comparable 
to a from Giglio et al. [2010] because here we used a subset 
of the active fires in each grid cell (only those associated 
with FG in ) to compute this ratio. We included the spatial 
buffer described above to allow for more conservative 
(smaller) estimates of FC out . Preliminary analysis suggested 
that increasing the buffer radius would reduce FC out and a, 
but would be at least partly compensated by increases in the 
7 term described below in our estimates of small fire burned 
area from equation (2). 

[ 16 ] To better understand the sensitivity of equation (2) to 
different satellite and environmental conditions, it can be 
rewritten in the following form, by replacing a with 
BA MCD(AAX IFC in and rearranging terms. 


where the total burned area, BA totah in each grid cell, i, month 
t, and aggregated vegetation class, v, is the sum of the original 
MCD64A1 burned area, BA MCD 64Al , and a new component 
associated with small fires, BA s f. As described below, our 
approach for estimating BA s f was statistical, and required 
information about the ratio of active fires in or near large 
burned areas to active fires outside of these areas. Our choice 
of the 0.25° grid and the averaging system for model para- 
meters described below reflected a balance between our goal 
of providing relatively high spatial resolution information for 
investigating regional to global scale Earth system science 
questions, and the requirements for maintaining a robust 
computation approach. 

[ 14 ] We estimated the small fire component as: 

BA sf (i, t, v) = FC ou ,(i , t, v) x a r , w x 7 ,. w , (2) 

where FC out is the total number of MOD14A1 active fires 
outside of burned areas in each 0.25° grid cell, a is the ratio 
of BA M cdm A \ to active fire counts within or near the perim- 
eter of these burned areas, and 7 is a scalar that modifies the 
burned area associated with FC out because these active fires 
often were associated with smaller bum scars compared to 
active fires within or near burned areas identified by the 
MODIS MCD64A1 product. Active fires within or near 
burned areas (the complement to FC out ) are hereafter referred 
to as FC in . a and 7 were separately estimated within each of 
14 continental-scale regions (r) for 3 seasonal periods ( 5 ), for 
5 aggregated vegetation classes (v), and for each year (v). 
A flow diagram showing how the different MODIS pro- 
ducts were used to compute equation (2) is provided in 
the auxiliary material (Figure SI). 1 

[ 15 ] We computed FC out by applying a 1-km buffer to the 
perimeter of MCD64A1 burned areas, and then counting the 
number of active fire detections outside of this area. FC in 
was defined as the complement to FC ot „, i.e., the active fires 
that overlap with MCD64A1 burned area or the surrounding 
buffer area. All the active fires within one month of the bum 
date were assigned to FC in . a had units of km 2 per active fire 
and was defined for our best case estimate using MCD64A1 
burned area and FC,„ derived from the Terra MOD14A1 
product. The spatial distribution of a was influenced by the 
efficacy of the burned area detection algorithm, the fre- 
quency of satellite overpasses, and aspects of fire behavior 


’Auxiliary materials are available in the HTML. doi:10.1029/ 
2012JG002128. 


BA s f(i, t, v) — BA MC d 64 ai x 


FC ou t(r,s,v,y ) 
FC in {r, s , v, y) 


x 7 j 


r,s,v{yi 


(3) 


Although this form was not used in our analysis and 
involves a different approach for averaging FC in and FC ou t, 
it is important for understanding the sensitivity of BA sf to 
different processes. Specifically, the general form of our 
approach depends on the ratio of active fires detected inside 
and outside of bum perimeters in each region and not the 
absolute number of active fires detected outside of bum 
scars. This may make the technique less sensitive to regional 
variations in satellite overpass dynamics and cloud cover, 
since the impacts of these factors would be expected to have 
similar effects for FCi„ and FC out . Hereafter, when we refer 
to FC out (and FC in ), we refer to the form of FC out ( i , t, v) 
described by equation (2). 

[ 17 ] The 7 in equation (2) was an additional unitless scalar 
applied to FC out to account for the observation that the set of 
active fires in FC out did not necessarily correspond to the 
same size distribution of bum scars as those in FCi„. Spe- 
cifically, in many regions FC out had smaller dNBR changes 
than FCin , indicating that on average, individual FC out active 
fires were associated with smaller bums. As described below 
in the results (section 3.5), FC out also had lower values of 
FRP, indicating that the active fires were smaller in size or 
cooler. 7 was estimated using the following equation: 

, , , _ dNBR out {r,s,v,y) - dNBR control (r,s,v,y) 

1 r,s,v ' y dNBR in {r,s,v,y) - dNBR control {r,s,v,y) ’ 


where dNBR ol „ was the difference normalized burn ratio 
corresponding to FC out , dNBR in was the same quantity but for 
FC in , and dNBR colltro! was derived using dNBR from nearby 
unbumed regions (for the same period, aggregated vegetation 
class, region, and year). We note that our approach assumed 
that dNBR differences between FC in and FC out , each relative 
to the control, were linearly proportional to subpixel burned 
area. While not strictly valid, this assumption was not entirely 
unreasonable provided the associated reflectance changes 
were not excessively large [ Roy and Landmann, 2005]. We 
also assumed that all of the dNBR differences were attribut- 
able to RH^rather than to a change in fire severity. This latter 
assumption was probably conservative with respect to the 
BA s f values we report in the results. Specifically, if bum 
severity was lower (per nr of burned area) for small fires 
outside of MCD64A1 burned areas, then the small fire burned 
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area fractions would have to be even higher to affect the same 
dNBR change (i.e., the 7 applied in equation (2) would need to 
be larger and thus closer to a value of 1). 

2.3. Global Application of Small Fire Burned 
Area Algorithm 

[is] We computed a and 7 for three different seasonal 
periods and five aggregated vegetation classes within 14 
continental-scale regions each year. We defined early, mid, 
and late season periods relative to the 2001-2010 mean peak 
fire month that was derived for each region using MOD14A1 
active fire observations. Early season was defined as 1— 
6 months prior to the peak fire month, middle season as the 
peak fire month, and late season as 1-5 months after the peak 
fire month. We chose this seasonal partitioning to separate 
active fires into three bins of roughly equal size, in terms of 
the number of active fires within each seasonal interval. 

[ 19 ] We also defined five aggregated vegetation classes 
within each continental region using the International Geo- 
sphere-Biosphere Programme (1GBP) vegetation classification 
scheme in the Collection 5 MCD12Q1 global land cover type 
product [Friedl et al. , 2002; Friedl et al . , 20 1 0] . All the forests 
(classes 1-5) were grouped together in one class as were 
closed and open shrublands, woody savannas, and savannas 
(classes 6-9). Grasslands (class 10) were considered another 
separate aggregated class. Croplands and cropland/natural 
vegetation mosaic types also were grouped together (classes 
12 and 14). The remaining vegetation types, including per- 
manent wetlands (class 1 1), urban or built up (13), snow and 
ice (15), and barren or sparsely vegetation (class 16) were 
combined into an ‘other’ aggregated vegetation class. We used 
2009 land cover type data from collection 5 for this mapping 
and anticipate improved performance with future versions of 
this product that remove known biases in high latitude regions. 

[ 20 ] The seasonal periods and aggregated vegetation clas- 
ses described above allowed a and 7 to capture some of the 
variability in environmental and ecosystem controls on fire 
behavior in each region, but were sufficiently coarse to allow 
for robust estimation of these parameters using available data. 

[ 21 ] We then computed BA s f using equation (2) with 
monthly varying values of FC out in each 0.25° grid cell and 
the matrices of a and 7 parameters described above that 
had dimensions of 10 years, 14 continental-scale regions, 
3 seasonal periods, and 5 aggregated vegetation classes. 

2.4. Estimation of Small Fire Burned Area During 
Gaps in the Burned Area (MCD64A1) Record 

[ 22 ] Application of equation (2) was possible only when both 
burned area tiles (MCD64A1) and active fires (MOD14A1) 
products were available. During 2001-2010, less than 1% of 
active fires occurred in places and months when monthly 
burned area tiles were unavailable. Thus, the approach 
described in the previous section (using equation (2)) was the 
most common approach for estimating BA s f during 2001- 
2010. Where MCD64A1 burned area tiles were not available 
(mostly in very high northern and southern latitudes and 
during a gap in June 2001), we could not estimate FC out , and 
thus BA s f in equation (2). For these periods and regions, we 
used the following extrapolation to estimate both BA sf and 
the total burned area. 

[ 23 ] First, we created matrices of the long-term mean values 
of a and 7 averaged using available observations during 


2001-2010. These matrices had dimensions of 14 continental- 
scale regions, 3 seasonal periods, and 5 aggregated vegetation 
classes. We hereafter refer to these climatological mean para- 
meters as a c and j c . We constructed an additional matrix, R out , 
of the mean ratio FC out to the sum of FC out and FC in with the 
same dimensions. To estimate burned area normally captured 
by MCD64A1 during months with gaps, we applied the fol- 
lowing equation: 

RTv/O/XU-I ] gap ( i- t, v) — F'C.'io/al ( f t, v) X (1 Rout (X A, v)) 

x a c (r,s,v), (5) 

where FC tota i is the total number of active fires observed in 
each 0.25° grid cell in a given month, separated by the five 
aggregated vegetation classes, 1 — R out is the mean fraction of 
FC tofa i that is within or near burned areas (this parameter 
varied by continental region, seasonal period, and aggregated 
vegetation class), and a c is defined in the preceding paragraph. 

[ 24 ] Similarly, we estimated the burned area associated 
with small fires during these gap periods using the following 
equation: 

BA s f gap (i, t, v) FCuoiaiU- t, v) x R nu t (V, s . vj 

x Oc(r,s,v) x J c (r,s,v), (6) 

where R ou , was the mean fraction of FC to tai associated with 
small fires (i.e., active fires outside of burned areas and vary- 
ing by continental region, seasonal period, and aggregated 
vegetation class) and a c and j c were climatological mean 
parameters that convert the product of FC, ota i and R out to 
burned area. Note that the product of FC to , a / and R oll , provides 
an estimate of FC ou , during gaps, so that equation (6) has 
a form similar to that of equation (2). This extrapolation 
approach accounted for less than 1% of the burned area in the 
2001-2010 S/fy time series presented below in section 3. 

2.5. Evaluation 

[ 25 ] We compared BA tota i from our approach with several 
regional burned area data sets derived from higher resolution 
satellite imagery and geographic information system (GIS) 
fire location and fire size data. In Alaska we used the point 
version of the Alaskan Large Fire Database (AKLFDB) 
developed by the Alaska Fire Service (http://fire.ak.blm.gov/ 
predsvcs/maps.php) [Kasischke et al., 2002] and the 
continental-scale polygon data set produced by the Monitoring 
Trends in Bum Severity (MTBS) project [Eidenshink et al., 
2007]. In Canada, we used the point version of the Canadian 
National Fire Database (CNFDB, http://cwfis.cfs.nrcan.gc.ca/ 
en_CA/nfdb/pnt, version 20101210) developed by the Cana- 
dian Forest Service [Amiro et al., 2001; Parisien et al., 2006]. 
We masked out Manitoba during 2001-2004, Nova Scotia 
during 2010, and Newfoundland and Labrador during 2001- 
2010 since fires from these provinces had not yet been added 
to the national database. For all of these comparisons, we first 
created raster maps of the observations on the same 0.25° grid 
as BA sf and BA lolal . 

[ 26 ] For the continental United States (CONUS), we sep- 
arated fires into several components. For wildland fires, we 
used observations from MTBS. To compare with our MODIS 
estimates, we regridded the National MTBS Burned Area 
Boundaries Data set (http://www.mtbs.gov/index.html) shape 
file perimeters to annual gridded maps. The MTBS protocol is 
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to map all fires on federal and state lands greater than 202 ha 
(500 acres) in the east and greater than 404 ha (1000 acres) in 
the west, with the dividing east-west line along the eastern 
border of North Dakota, South Dakota, Nebraska, Kansas, 
Oklahoma, and Texas. Although MTBS does not exclude 
prescribed fires with this protocol, these contributions are rel- 
atively small and thus cumulative MTBS burned area estimates 
for the United States are very similar to wildland fire estimates 
reported by National Interagency Fire Center (MFC) 
[Eidenshink et al., 2007]. For agricultural burned area, we 
compared our estimates with the cropland product developed 
by McCarty et al. [2009]. The McCarty et al. [2009] estimates 
were available during 2003-2007 and were derived using 
regionally specific dNBR thresholds that were optimized based 
on 296 ground observations (GPS bum perimeters) from 
burned and plowed fields. To compare with these observations, 
we first calculated the amount of BA totaI , BA s f, and BA mcd64ai 
in CONUS that occurred in croplands and cropland/natural 
vegetation mosaic types using the 500 m MCD12Q1 global 
land cover product. 

[27] In a final step of evaluation in CONUS, we compared 
our estimates with prescribed fires from NIFC [NIFC, 2012], 
For this, we computed a residual component after removing 
burned area overlapping with MTBS fire scars and burned 
area in agricultural regions, following the methods described 
above. The NIFC prescribed fire statistics were available for 
individual states, and we compared these estimates with the 
MODIS-based burned area for CONUS as a whole and on a 
state-by-state basis. 

[28] We also compared our estimates of BA MCD64Al and 
BA tota i with aggregated burned area statistics derived from 
Landsat imagery from the National Forest service in Portugal 
(www.afh.min-agricultura.pt/portal/dudf/cartografia/cartograf- 
areas-ardidas- 1990-2009) [e.g., Marques et al., 2011] and a 
savanna region in southern Mali, Africa based on the study of 
Laris [2005], For the Portugal data set, we created raster maps 
of the polygon information on the same 0.25° as BA s f. For the 
Mali comparison, we extracted and averaged together all of 
the 0.25° grid cells during the 2002-2003 fire season within 
the Landsat ETM + scene WRS- 199/52 that was the focus 
of the Laris [2005] analysis. In our validation, we attempted to 
use all publicly available regional burned area data that we 
could find. We recognize that these data sets are biased toward 
North America and anticipate increased availability of high 
quality burned area observations from many other continents in 
the future from Landsat as a part of ongoing Global Observa- 
tion of Forest and Land Cover Dynamics (GOFC-GOLD) 
efforts and other activities. 

2.6. Sensitivity Tests 

[29] We also separately compared MODIS fire radiative 
power (FRP) of active fires outside of burned areas ( FC out ) 
with those within or near burned areas ( FC in ). Fire radiative 
power provides information about some combination of fire 
temperatures, flaming versus smoldering fire stages, and fire 
area within a 1-km MODIS pixel [ Freeborn et al, 2008; 
Ichoku et al., 2003; Wooster et al, 2003]. As described 
below in the results, these FRP comparisons indicated that 
active fires outside of burned areas were smaller and/or 
cooler. We did not use FRP directly to estimate 7 because it 
remains uncertain how instantaneous fire behavior may be 
related to the final size of an individual fire. 


[30] We also conducted several tests to evaluate the sen- 
sitivity of our results to different assumptions with our 
approach. These included 1 ) using 8-day surface reflectance 
(MOD09A1) to estimate dNBR instead of the 16-day vege- 
tation index product described above; 2) constructing the 
dNBR control using a temporal sampling approach instead 
of the spatial sampling approach described above, 3) repla- 
cing Terra with Aqua active fires for the estimation of FC in , 
FC out , and a, 4) the same as 3 but also using Aqua surface 
reflectance to compute 7 ; and 5) only using nominal and 
high confidence active fire data to determine FC out , a, and 7 . 
For the temporal sampling estimate of dNBR co „ lro/ , we 
obtained this estimate by sampling within burned area peri- 
meters in the year prior to burning, and also by stratifying by 
aggregated vegetation class and seasonal period. For all five 
of these tests, we compared burned area estimates with 
results from our primary approach for the year 2010. 

2.7. Carbon Emissions Modeling 

[ 31 ] To test the impact of adding small fire burned area on 
emissions, we ran the GFED3 biogeochemical model [van 
der Werf et al., 2010] with BA MCD64Al as well as with 
BA tota! . Boosting burned area will increase emissions, but the 
relation is not necessarily linear for a number of reasons. 
First, it depends on where and when the small fires occur 
because 1) fuel loads vary considerably across the landscape 
as a function of vegetation type and 2) the fraction of fuels 
consumed (combustion completeness) varies through time as 
a function of fuel moisture levels. Second, if the increase in 
burned area occurs in areas with frequent fires, then the fuel 
loads may not have had time to recover from the previous 
fire, leading to proportionally smaller emissions increases. 

[32] We completed two simulations at a 0.25° spatial reso- 
lution. All input data sets except burned area, fire persistence, 
and the fraction of burned area occurring in tropical peatlands 
were the same and were derived from the original 0.5° GFED3 
runs described by van der Werf et al. [2010], Within each 
0.25° grid cell, we kept track of the sub-grid partitioning of the 
500 m burned area for BA MCD 64A1 over different land cover 
classes, tree cover bins, and peatland maps [van der Werf et al. , 
2010]. We did the same with BA s f by tracking the vegetation 
types corresponding to the active fires in FC out , weighted by 
the appropriate a and 7 which were unique for different 
aggregated vegetation classes, regions, and years. 

[33] We made several small changes to the biogeochemi- 
cal model compared to van der Werf et al. [2010]. Here we 
used 2001-2010 burned area observations to spin up the 
biogeochemical model. Also, in the original version of 
GFED3, burned area in forests in deforestation regions was 
boosted with fire persistence to gain consistency with pub- 
lished deforestation assessments [ERODES, 2012]. Here we 
followed the same approach and applied an additional linear 
scalar of 1.42 for the simulation with BA MCD64Al derived 
from a similar comparison as in van der Werf et al. [2010] 
but now with more years of deforestation data in the Brazi- 
lian Amazon from the Amazon Deforestation Monitoring 
Project [ ERODES, 2012] in the optimization. For the BA tota i 
run, we did not apply this scale factor to burned area in 
deforestation regions as deforestation rates predicted using 
this approach were larger than several recent estimates 
[PRODES, 2012; Hansen et al, 2008]. Instead, when burned 
area exceeded the observed deforested area we capped the 
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Figure 1. The number of active fires within or near burned areas ( FC in , red bars) or outside of burned 
areas ( FC out , light orange bars) during 2001-2010. The burned area observations used to define FC in 
and FC out are from the 500 m MCD64A1 product, and data on the active fires are from the Terra 
MOD14A1 product. The FC out estimates are stacked on top of the FC in estimates, so their sum (the height 
of the total bar) is equal to FC tota i. Different continental-scale regions are abbreviated as follows: Boreal 
North America (BONA), Temperate North America (TENA), Central America (CEAM), Northern Hemi- 
sphere South America (NHSA), Southern Hemisphere South America (SHSA), Europe (EURO), Middle 
East (MIDE), Northern Hemisphere Africa (NHAF), Southern Hemisphere Africa (SHAF), Boreal Asia 
(BOAS), Central Asia (CEAS), Southeast Asia (SEAS), Equatorial Asia (EQAS), and Australia (AUST). 
A map of these regions is given in van der Werf et al. [2010]. 


burned area associated with deforestation at similar levels as 
in the BA mcdmax run and applied this amount to the defor- 
estation module in GFED3. The excess burned area in 
tropical forests was treated separately as a class of fires not 
directly associated with land clearing (i.e., as escaped fires 
or forest degradation fires). 

[ 34 ] We note that even though small fires are probably less 
severe, as described below in Section 3 in terms of dNBR 
and FRP observations, this potential bias was taken into 


account in our emissions modeling approach because BA s f 
was computed using 7, which in turn depends on dNBR. 

3. Results 

3.1. Spatial Distribution of Active Fires 
and Model Parameters 

[ 35 ] We found good agreement between the spatial dis- 
tribution of active fires and burned area in boreal forests of 
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Figure 2. Annual active fire counts from MODIS Terra (number of fire counts per 0.25° grid cell per 
year), (a) Active fires within the perimeter of the MODIS burned area product MCD64A1 (same as FC in 
in the text), (b) Active fires outside of MODIS burned area (same as FC out ). Figure 2b is the complement 
to FC in in Figure 2a. (c) The total number of active fires in each grid cell ( FC tota i ). (d) the fraction of active 
fires outside of bum perimeters to total active fires (FC out /FC tota/ ). All of the panels show annual means 
constructed using 2001-2010 observations. 


North America and in savannas across Africa and Australia. 
In these regions most active fires were within or near bum 
scars, and so FC in was larger than FC ou , (Figures 1 and 2). 
In contrast FC ou , was considerably larger than FC in in tem- 
perate North America, Europe, Central America, South 
America, Southeast Asia, and equatorial Asia. As a function 
of vegetation type, on all continents FC out was larger than 
FC in in croplands and in cropland/natural vegetation mosa- 
ics. This was expected given the small field sizes and frag- 
mented landscape that exist in many agricultural regions 
[. McCarty et al., 2009] and the challenge of detecting bum 
scars in areas undergoing other forms of land management 
such as harvesting, plowing, or irrigation. In many tropical 
regions FC out was a larger fraction of the total number of 
active fires at the end of the fire season (Figure 1). 

[ 36 ] a, the ratio of burned area to active fires, varied 
considerably as a function of vegetation type. In boreal for- 
ests, a was often less than 1 km 2 /active fire, probably as a 
result of multiple satellite overpasses at higher latitudes each 
day and high fuel loads leading to multiple active fire 
detections for a given burned area perimeter (Figure 3). 
In many tropical savanna regions, in contrast, a often varied 
between 2 and 5 km 2 /active fire, probably a result of faster 
moving (and short-lived) grass fires, relatively low fuel 
levels, and longer gaps between successive overpasses that 
limited the number active fires observed within bum peri- 
meters [ Giglio et al., 2010], 

[ 37 ] 7, the scalar adjustment to a for active fires outside of 
burned areas, also varied considerably as a function of 


region and vegetation type (Figure 3). This scalar typically 
varied between 0.2 and 0.6 for different seasonal periods and 
vegetation types, with larger values observed for Northern 
Flemisphere Africa, Central Asia, Southeast Asia, and 
smaller values in Temperate North America and Europe. 
Example distributions of normalized bum ratio used to 
derive 7 for the most commonly burning aggregated vege- 
tation class in each region, during the peak month of burn- 
ing, are shown in Figure 4. 

[ 38 ] In terms of spatial variability of BA sf across the dif- 
ferent continental-scale regions shown, regional differences 
in FC out contributed to 43% of the variance, a contributed 
36%, and 7 contributed to 21%. Considering the full range 
of mean model parameters, FC out varied by a factor of 26 
across the different continental-scale regions, FC oll /FC tota! 
varied by factor of 3.8, a varied by a factor of 5.0, and 7 
varied by a factor of 2.8. 

3.2. Comparison With Regional Burned Area Products 

[ 39 ] In Alaska, including small fires increased the MODIS 
estimates of burned area by about 8% during 2002-2010, 
bringing it into closer agreement with AKLFDB and MTBS 
products (Figure 5a). The original BA MCD 64AI had a mean 
burned area of 0.65 Mha/yr during 2002-2010, approxi- 
mately 24% lower than the mean from AKLFDB. The small 
fire algorithm (BA s y) added about 0.06 Mha/yr, so that the 
total burned area from MODIS was 17% lower than 
AKLFDB. In all years, the spatial correlation coefficient 
between BA MCD 64A1 and AKLFDB increased by a small 
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Figure 3. Model parameters a and 7 for different continental-scale regions and aggregated vegetation 
classes averaged over 2001 - 2010 . 


amount when the small fire component was added. Averaged 
over each year during 2002 - 2010 , the mean spatial correla- 
tion coefficient between BA mcd64a 1 and the AKLFDB was 
0.69, and this increased to 0.71 for BA tota i. For this same 
period, fires less than 1000 ha in size accounted for 1.5% of 
burned area in the AKLFDB. For this size class (i.e., for grid 
cells that had more than 50% of total burned area derived 
from this size class), the BA s f algorithm increased MOD1S 
burned area by 42%, whereas for grid cell dominated by large 
fires greater than 1000 ha, BA s f increased MODIS burned 
area by only 8 %. Further visual inspection of the spatial 
pattern of the FC out set of active fires in Alaska showed that 
almost all were within or near (less than 1 km) of AKLFDB 
polygons, indicating that most of these pixels were associated 
with existing fires in the database, but in areas that may have 
burned only partially, as expected for example for areas near 
the fire perimeter or near interior unbumed islands. We note 
that in comparisons with geographic information system 
polygon products, 500 m burned area products such as 
MCD64A1 are expected to yield lower estimates because 
small unbumed islands will be excluded from this product 
when the algorithm is applied on a pixel by pixel basis. For 
example, using a dNBR-based classification tailored specif- 
ically for fires in Alaska, Beck et al. [2011] estimated that 
unbumed islands accounted 3% of the interior area of low 
severity polygons and 7% for high severity polygons. 


[ 40 ] In Canada, the CNFDB was higher than the original 
BA MC d 64 A \ in all years. The small fire algorithm increased the 
MODIS mean during 2002-2010 (Figure 5b) and slightly 
increased the spatial correlation in all but one of these years in 
non-agricultural areas. BA S/ considerably modified MODIS 
burned area estimates across the wheat belt in Saskatchewan 
and across the southwestern parts of British Columbia and 
Alberta (Figure S2 in the auxiliary material). Burned area more 
than tripled in agricultural areas, from 0.04 to 0.14 Mha/yr, 
when small fires were included. In Canada, fires less than 
1 000 ha in size accounted for 4. 1 % of burned area in the point 
version of the CNFDB. For this size class, the BA s j algorithm 
increased MODIS burned area by 98%, whereas for larger 
fires greater than 1000 ha, BA s f . increased MODIS burned area 
by only 4%. 

[ 41 ] In the continental U.S. the small fire algorithm 
increased the total amount of MODIS burned area consid- 
erably, by approximately 75% during 2001-2010 (Table 1), 
and brought the MODIS-based estimates into closer agree- 
ment with several estimates. For areas that experienced large 
wildland fires, defined as 0.25° grid cells that had at least 1 
fire during 2002-2010 in the MTBS database, the small fire 
algorithm increased MODIS burned area by 16% (Figure 5c). 
BA MC d 64 ai was 60% of the MTBS mean during 2002 - 2010 , 
and this percentage increased to 70% for BA tota/ . For the 
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Figure 4. The distribution of difference normalized bum ratio (dNBR) for active fires within burned 
areas from MCD64A1 (red), outside of burned areas (orange) and for control areas (blue). The distribu- 
tions were constructed during the peak fire month using the most commonly burning vegetation type with 
each region. The distributions were generated using observations made between 2001 and 2010. 


MTBS data set, fires less than 1000 ha in size accounted for 
12% of the total burned area. For this size class, the BA s f 
algorithm increased MODIS burned area by 70%, whereas 
for larger fires greater than 1000 ha, BA s f increased MODIS 
burned area by only 12%. 

[ 42 ] Both of the MODIS estimates were broadly consistent 
with MTBS burned area given, for example, that 28% of the 
area within MTBS polygon perimeters was classified as 
“unbumed to low severity” for the Pacific Northwest and 
Pacific Southwest during 1984-2005 [Schwind et al., 2008], 
It was not possible to further quantify if BA s f improved 
MODIS estimates for MTBS because statistics on unbumed 
islands were not readily available for the continental-scale 
MTBS polygon data set, and because it was not possible to 
simulate the “unbumed to low severity” MTBS class with 
our BA MC d64a\ and BA s f algorithms, given that these classes 
were computed from different temporal intervals of dNBR 
for each fire in the MTBS database. Also, it is likely that 
some fires not digitized as a part of the MTBS project were 
included in some of the overlapping 0.25° grid cells for 
BAtotal • 


[ 43 ] For croplands in CONUS, including small fires 
improved comparisons with published estimates for this 
biome. The original BA MCD 64A1 had a 2003-2007 mean of 
0.26 ± 0.06 Mha/yr for U.S. croplands, defined by burned 
area in cropland and cropland and natural vegetation mosaic 
500 m pixels using the IGBP classification from Friedl et al. 
[2010], The small fire algorithm tripled this estimate, 
yielding a total burned area mean of 0.85 ± 0.08 Mha/yr and 
bringing our estimates into closer agreement with previous 
estimates for this biome of 1.24 ± 0.07 Mha/yr from 
McCarty et al. [2009] for this same period (Figure 5d). 

[ 44 ] After removing the wildland fire and agricultural 
burned area components described above from BA M cdma\ 
and BA tota i, we obtained residuals for CONUS of 0.22 ± 
0.06 Mha/yr Mha/yr and 0.69 ± 0.28 Mha/yr during 2002- 
2010. Comparing these residuals to the NIFC reports of 
prescribed fires for this period (1.03 ± 0.15 Mha/yr), small 
fire burned area considerably improved the magnitude of the 
agreement though BA tota i remained below NIFC prescribed 
fire reports for most years (Figure 5e). The spatial distribu- 
tion of this residual burned area for BA total was concentrated 
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Year 


Figure 5. Comparison of MODIS burned area products 
with regional burned area products for North America. The 
regional observations are shown in black, the MODIS 
MCD64A1 burned area is shown in red, and MODIS total 
burned area, including the small fire burned area developed 
here, is shown in blue. In Figure 5a, observations for Alaska 
are provided from both the Alaska Large Fire Database 
(AKLFDB, black squares) and the Monitoring Trends in 
Bum Severity (MTBS) project (black circles). 

in the Southeast (Figure 6), which is consistent with exten- 
sive use of prescribed burning for forest and plantation 
management in this region. Small fires modestly improved 
the spatial agreement across the 48 states between the NIFC 
prescribed fire estimates and the MODIS residual, with the 


correlation coefficient for the 2002-2010 mean increasing 
from 0.28 for BA MCD64Al to 0.30 for BA total . 

[ 45 ] For Portugal, for the entire domain of the country the 
two MODIS burned area estimates were lower than the 
national polygon data set during 2002-2009, by 24% for 
BA MCD64A i and 7% for BA tota/ (Figure S3 in the auxiliary 
material). In 7 of 8 years, BA tota i had a higher spatial corre- 
lation with the national data set than BA MC dma\, an d over 
the full period small fires increased the spatial correlation 
coefficient on average from 0.67 to 0.70 (Figures S4 and S5 
in the auxiliary material). 

[46] For the savanna ecosystem region in Mali studied by 
Laris [2005], including small fires increased MODIS burned 
area during the 2002-2003 fire season from 24% to 33% of 
the study area. Compared to the Landsat-derived estimate of 
57% for this domain [Laris, 2005], accounting for small fires 
improved the agreement but did not eliminate the low biases 
observed for the MODIS burned area product. The relative 
monthly increases from including small fires were largest at 
the beginning and end of the fire season, when the preva- 
lence of small and fragmented bums was higher based on the 
Landsat-derived perimeter to burned area ratios (Table 2). 

3.3. Impacts of Small Fires on Global and Regional 
Burned Area Patterns 

[ 47 ] At a global scale, including small fires increased 
burned area by 35% or 120 Mha/yr during 2001-2010 
(Table 1). Most of the global increase was driven by small fires 
in Southern Flemisphere Africa (42.3 Mha/yr), Northern 
Hemisphere Africa (36.0 Mha/yr), Southern Hemisphere 
South America (10.2 Mha/yr), Central Asia (10.1 Mha/yr), 
and Southeast Asia (6.0 Mha/yr). In terms of relative changes, 
Equatorial Asia had the largest increase (157%), followed by 
Central America (143%), Europe (112%), Southeast Asia 
(90%), and the Middle East (90%). As a function of vegetation 
type, woody savannas (46.4 Mha/yr), savannas (30.5 Mha/yr), 
croplands (11.8 Mha/yr), cropland and natural vegetation 
mosaics (1 1 .2 Mha/yr), and grasslands (9.6 Mha/yr) accounted 
for most of the increase in burned area (Table 3). Averaged 
over the decade we estimated about 464 Mha, or 3.5% of the 
ice-free land surface, burned each year. 

[48] At a regional scale, large relative increases in burned 
area occurred across eastern North America, southern Mex- 
ico, Central America, Eastern Europe, Central Asia, and 
many countries in Southeast Asia (Figures 7 and 8). Many of 
these increases were driven by a much higher detection 
efficacy of the small fire burned area algorithm in croplands, 
cropland and natural vegetation mosaics, and evergreen 
broadleaf forests (Table 3). From an emissions perspective, 
the increase of 119% of burned area in evergreen broadleaf 
forests is particularly important because of the high fuel 
loads present in this biome. Equatorial Asia had one of the 
largest relative increases in burned area of any region 
(157%). The small fire burned area product in this region 
revealed a much more extensive distribution of burning 
across Sumatra and in western and interior areas of Borneo 
(Figure 8). Although the amount of burned area in equatorial 
Asia is quite low compared to Africa and other regions, 
many of these fires occur in fuel-rich tropical peatlands 
[Page et al., 2002] that have a high sensitivity to ENSO- 
induced variability in regional climate [Field et ah, 2009; 
Tosca et al., 2011]. 
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Table 1. Burned Areas by Region Averaged During 2001-2010 


Region 3 

Original MODIS 
500 m 

Burned Area 
( BAmcd64ai\ Mha/yr) 

Burned Area 
From Small Fires 
(. BA s f, Mha/yr) 

Total Burned Area 
(BA lota ,; Mha/yr) 

Increase From 

BA sf (%) 

CV b for 

BAmCD64A1 

CV for 
BA s f 

CV for 
BAtotaI 

BONA 

2.3 

0.3 

2.6 

14 

0.56 

0.35 

0.50 

TENA 

1.6 

1.2 

2.7 

75 

0.41 

0.35 

0.30 

CEAM 

1.4 

2.0 

3.4 

143 

0.51 

0.23 

0.30 

NHSA 

2.9 

2.3 

5.1 

79 

0.22 

0.18 

0.19 

SHSA 

18.5 

10.2 

28.7 

55 

0.47 

0.20 

0.36 

EURO 

0.5 

0.6 

1.1 

112 

0.42 

0.33 

0.31 

MIDE 

0.7 

0.6 

1.2 

90 

0.40 

0.37 

0.38 

NHAF 

117.9 

36.0 

154.0 

31 

0.10 

0.11 

0.10 

SHAF 

126.2 

42.3 

168.5 

34 

0.07 

0.05 

0.05 

BOAS 

6.0 

3.8 

9.8 

62 

0.76 

0.47 

0.59 

CEAS 

14.1 

10.1 

24.2 

72 

0.31 

0.17 

0.22 

SEAS 

6.7 

6.0 

12.6 

90 

0.37 

0.28 

0.30 

EQAS 

0.7 

1.1 

1.8 

157 

0.82 

0.42 

0.52 

AUST 

45.2 

3.3 

48.5 

7 

0.49 

0.42 

0.49 

Global: 

344.6 

119.7 

464.3 

35 

0.07 

0.04 

0.06 


“The different regions are defined in the caption for Figure 1. 

b CV is the coefficient of variation computed over the 2001 2010 period using the annual mean observations from each year. CV is defined as the 
standard deviation divided by the mean. 


3.4. Small Fires Effects on Seasonal and Interannual 
Variability of Burned Area 

[ 49 ] In tropical forest and savanna regions, including small 
fires modified the annual cycle of burned area by extending 
the fire season. Specifically, small fires increased burned 
area by a larger amount at the end of the fire season than at 
the beginning of the fire season (Figure 9). This effect was 
most pronounced in Central America, Equatorial Asia, 
Southeast Asia, and South America. The greater contribution 


from active fire detections to BA totaI in the late dry season 
may also reflect difficulties with burned area mapping dur- 
ing the dry-wet transition due to higher cloud cover or 
aerosol concentrations. In Central and Boreal Asia, in con- 
trast, small fires had the largest impact on total burned area 
during spring, at the onset of the fire season. These changes 
in the annual cycle were mostly driven by patterns of FC out 
and are visible in the ratio of FC out to FC tota i shown in 
Figure 1. Seasonal changes in a and 7 appeared to have only 



(a) Total burned area 


(b) Wildland fire burned area (burned area in MTBS grid cells) 

-120 -110 -100 -90 -80 -70 



(c) Agricultural burned area 



(d) Prescribed fire and residual burned area 


0.01 


0.1 


1 10 100%/yr 


Figure 6. MODIS burned area fraction for the continental U.S. (. BA lola i) for (a) all fires, (b) burned area 
annually coincident at a 0.25° spatial resolution with Monitoring Trends in Bum Severity (MTBS) project fire 
scars, (c) agricultural fires, and (d) prescribed/other fires. Figure 6c was constructed by sampling BA tota / in 
areas where 500 m vegetation cover maps indicated either croplands or croplands and natural vegetation 
mosaics. Figure 6d was constructed by subtracting fires shown in Figures 6b and 6c from those shown 
in Figure 6a. This residual includes prescribed fires and other fires that were not associated with MTBS 
perimeters or cropland areas. 
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Table 2. MODIS Burned Area Comparisons With Landsat-Derived Estimates From Laris [2005] 


Burned Area Estimates for Southern Mali, 

MODIS-Derived Burned Area Africa [From Laris , 2005] 


Time Period 

BA MC d6 4A i 
(% of Study Area 
per Month or Year) 

BA tota i 

(% of Study Area 
per Month or Year) 

Increase From 
Small Fires (%) 

Landsat-Derived Burned Area 
(% of Study Area 
per Month or Year) 

Landsat-Derived Perimeter 
of Burned Area Ratio 

September-November 2002 

i.i 

4.9 

345 

17.6 

0.305 

December 2002 

14.3 

17.9 

25 

23.9 

0.235 

January 2002 

6.6 

7.5 

14 

10.5 

0.217 

February 2002 

2.0 

2.6 

30 

4.9 

0.429 

Total 2002-2003 fire season 

24.0 

32.9 

37 

56.9 



a secondary effect on the annual cycle of small fire burned 
area, and also on patterns of interannual variability. 

[so] One possible mechanism for the larger relative con- 
tribution of small fires to total burned area at the beginning 
and end of the fire season is that fires during these periods 
are more fragmented as a consequence of higher vegetation 
and soil moisture levels. During these times, patches of the 
landscape may have fuel moisture levels that are too high to 
bum, limiting the spread of large fires. This hypothesis is 
consistent with the mechanisms proposed by Laris [2005] to 
describe the seasonal dynamics of fires in wooded savanna 
ecosystems in Mali. Fire perimeter to burned area ratio 
observations from Laris [2005] indicate that larger fires 
occur in the middle of the fire season, when fuel moisture 
levels are sufficiently low to allow combustion in all vege- 
tation types present in landscape mosaics (Table 2). 

[ 51 ] From year to year, small fire burned area was con- 
siderably less variable than MCD64A1 burned area in many 
regions, including Central America, Southern Hemisphere 
South America, Southeast Asia, and Central Asia (Figure 10 
and Table 1). Globally, the coefficient of variation for BA s f 
was approximately 40% lower than that for BA MCD6 4A1 
(Table 1). One hypothesis for the less variable BA s f is that 
these fires are more closely linked with different forms of 
land management, such as agricultural waste burning or 
prescribed burning in forest plantations. Specifically, 
incentives may exist for land managers not to allow these 
fires bum out of control and damage nearby areas. If and 
when these fires do escape control, it may be as a conse- 
quence of anomalously dry climate conditions, leading to 


larger and more severe fires, and ones that would be more 
easily detected by the BA MCD 64AI product. 

3.5. Sensitivity Test Results 

[ 52 ] Mean fire radiative power for individual fire detec- 
tions was lower on average for active fires outside of 
MCD64A1 burned area than for active fires within the 
perimeter of these bum scars (Table 4). This is consistent 
with the premise that active fires associated with FC out were 
often smaller and thus more difficult to detect using burned 
area algorithms applied to moderate resolution surface 
reflectance time series. In many regions, the cumulative fire 
radiative power outside of burned areas, defined as the mean 
fire radiative power per active fire detection multiplied by 
the total number of active fires in this set during the time of 
the Terra overpass, was considerably larger than the cumu- 
lative fire radiative power within burned areas. 

[ 53 ] Sensitivity simulations for the year 2010 provided 
some information about the likely uncertainty in our esti- 
mates of burned area associated with small fires. Using 
active fires from Aqua increased our small fire burned area 
estimates by approximately 49% at a global scale. Two 
thirds of this increase was driven by a greater fraction of 
active fires occurring outside of burned areas during Aqua 
overpasses than during Terra overpasses. A likely explana- 
tion for this is that the diurnal cycle for smaller fires has a 
stronger peak in mid-day and thus proportionally more are 
detected during the nominal 1:30 pm overpass for Aqua. 
This finding is also consistent with diurnal cycles derived for 
different vegetation types from GOES that show more pro- 
nounced diurnal cycles in agricultural regions and tropical 


Table 3. Burned Area by Vegetation Type Averaged During 2001-2010 


Vegetation Type 

Original MODIS 500 m 
Burned Area 
( BA MCD 64Ai\ Mha/Year) 

Burned Area From 
Small Fires 
{BA s j\ Mha/Y ear) 

Total Burned Area 
( BA totab Mha/Year) 

Increase From 
Small Fires (%) 

Evergreen Needleleaf Forests 

1.7 

0.3 

2.0 

17 

Evergreen Broadleaf Forests 

3.1 

3.6 

6.7 

119 

Deciduous Needleleaf Forests 

1.1 

0.2 

1.3 

22 

Deciduous Broadleaf Forests 

3.4 

1.3 

4.7 

39 

Mixed Forests 

1.1 

0.8 

1.9 

75 

Closed shublands 

2.1 

0.7 

2.8 

33 

Open shrub lands 

17.9 

1.7 

20.0 

10 

Woody Savannas 

136.2 

46.4 

182.6 

34 

Savannas 

114.3 

30.5 

144.8 

27 

Grasslands 

38.5 

9.6 

48.1 

25 

Permanent Wetlands 

1.2 

0.8 

2.0 

67 

Croplands 

9.6 

11.8 

21.4 

123 

Urban 

0.1 

0.5 

0.5 

834 

Croplands and Natural Vegetation Mosaics 

14.3 

11.2 

25.5 

79 

Barren 

0.2 

0.1 

0.3 

71 
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Figure 7. Percentage per year of burned area fraction in each 0.25° grid cell associated with (a) the orig- 
inal 500 m MODIS burned area product (BA mcd(AA \), (b) the small fire burned area product developed 
here ( BA s j ), and (c) the total MODIS burned area ( BA tota i )■ (d) The ratio of small fire to total burned area 
in each grid cell ( BA s f/BA tota i ). These panels show the annual mean burned area fractions averaged using 
all available data during 2001 - 2010 . 


forest ecosystems than for fires occurring in temperate and 
boreal shrublands and forests [Mu et al., 201 1], The inferred 
values of 7 also increased, suggesting the size of individual 
fire scars for FC out derived from Aqua also increased. 

[54] Sampling bum perimeters in the exact same location but 
in the year before burning as a dNBR control in equation (4) 
reduced our estimates of small fires by approximately 15% 
globally. Similarly, using 8 -day surface reflectance (MOD09A1) 
instead of the 16-day vegetation index product (MOD13A1) 
to derived dNBR reduced our estimates by 16% globally. One 
possible issue related to use of the 8 -day surface reflectance is 
the presence of residual cloud contamination [Loboda et al., 
2007] that may be minimized with the compositing and lon- 
ger sampling interval of the MOD13A1 product. Using only 
middle and high confidence active fire retrievals from Terra 
for the computation of a and 7 had only a minor effect on our 
results, increasing small fire burned area by 6 %. Within indi- 
vidual regions, the impact of changing these model parameters 
varied considerably (Table SI in the auxiliary material). 
An important next analysis step is to assess how the different 
assumptions examined in the sensitivity analysis impact the 
comparisons with evaluation data sets from the Canada, the 
United States, Portugal, and Mali. This may guide future 
refinements to estimates of small fire burned area. A summary 
of the sensitivity results and their impact on BA tota i for differ- 
ent regions (and globally) is shown in Table SI. 

[55] Among the different combinations of sensitivity tests, 
using middle and high confidence active fires and observa- 
tions from Aqua would lead to a 55% increase in global BA s f 
(assuming the changes from these two tests are additive) 


whereas using Terra active fires with 8 -day surface reflec- 
tance and an alternative approach for constructing the pre- 
fire dNBR control would lead to a 31% decrease. This —31 
to +55% range across all the tests provided some informa- 
tion about uncertainty levels. More formal estimates of 
uncertainties (i.e., 95% confidence limits) were not possible 
at this stage for BA sf because of a lack of quantitative 
information regarding uncertainties in the different satellite 
products we used and also because of a paucity of high 
resolution burned area observations required for validation. 
Approaches for quantifying and reducing uncertainties in 
future work are discussed in Section 4. 

3.6. Effects of Small Fires on Biomass 
Burning Carbon Emissions 

[56] Global carbon emissions during 2001-2010 increased 
by 35%, from 1.9 Pg C/yr to 2.5 Pg C/yr, when BA s/ was 
included in GFED3 simulations (Figure 11 and Table 5). 
Within most regions, the relative increase in emissions from 
including small fires was smaller than that for burned area, 
as a consequence of more BA ^occurring in vegetation types 
with lower fuel loads (e.g., agriculture) and during times 
when climate and fire persistence controls on combustion 
completeness caused lower levels of fuel consumption 
within the model. One exception to this was Central Asia 
where small fire burned area increased carbon emissions 
from vegetation types with higher fuel loads. Even though 
13 of 14 regions had relative increases in emissions that 
were smaller than the relative increases in burned area, 
globally the percent increase in emissions (35%) was the 
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MCD64A1 burned area 


Small fire burned area 


(a) 




(c) 


0.01 0.1 1 10 100%/year 

Figure 8. (left) The MOD1S burned area product ( BA mcd64ai ) and (right) the small fire product devel- 
oped here ( BA s j ) for temperate North America, Central America, and Indonesia. Units are percentage of 
each 0.25° grid cell per year averaged over 2001-2010. 



same as for burned area. This was possible because large 
contributions to the total carbon emissions from BA s f origi- 
nated from regions with low levels of burned area, but high 
levels of mean annual emissions. Specifically, Central 
America, Northern Hemisphere South America, Southern 
Hemisphere South America, Southeast Asia, and Equatorial 
Asia accounted for over 42% of the carbon flux derived from 
small fires, yet only 1 8% of BA s f. The standard deviation of 
annual mean carbon emissions from these regions during 
2001-2010 increased as a consequence of small fires 
increasing burned area by a considerable amount in tropical 
forest areas. In other regions, with the exception of Northern 
Hemisphere Africa, the standard deviation of annual carbon 
emissions also increased during 2001-2010, though by a 
smaller amount. 

4. Discussion 

4.1. Overview 

[ 57 ] The difficulty of detecting small fires with moderate 
resolution surface reflectance imagery has been noted in 


many studies \Eva and Lambin, 1998; Simon et al., 2004; 
Hoelzemann et al., 2004; Laris, 2005; Roy and Landmann, 
2005; Silva et al., 2005; Loboda et al., 2007; McCarty 
et al., 2009; Roy and Boschetti, 2009]. Burned area detec- 
tion algorithms have to distinguish fires from many other 
forms of land cover change, including seasonal changes in 
leaf area and leaf senescence, synoptic and seasonal changes 
in soil moisture, and other forms of disturbance such as 
harvesting, plowing, insect outbreaks, and hurricanes. This 
classification problem places stringent requirements on the 
thresholds for change detection. As a consequence, a sub- 
stantial fraction of a pixel must be completely burned to 
avoid serious commission errors that would degrade the 
usefulness of a regional or global burned area product 
derived using moderate resolution surface reflectance 
imagery. Indeed, many developers of global burned area 
products have opted to accept larger omission errors than 
commission errors in their algorithm design [Roy and 
Boschetti, 2009]. We argue that while it will be possible to 
refine burned area detection algorithms in future work, the 
intrinsic tradeoff between commission errors and threshold 
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Figure 9. Seasonal dynamics of burned area for different continental scale regions (with units of Mha per 
month) for the original 500 m MODIS burned area product ( BA mcdma i ) denoted by red bars and the addi- 
tional small fire burned area product developed here (BA s j) denoted by orange bars. These panels show the 
annual mean burned area averaged during 2001-2010. 


relaxation may make it impossible for these algorithms to 
work effectively for some sub-pixel bums. Instead, other 
approaches are needed to quantify burned area and emis- 
sions for this class of fires. 

[58] Here we developed an approach for quantifying these 
fires consisting of two stages: (1) we computed number of 
active fires inside and outside of burned areas (FC in and FC out ) 
to identify the relative contribution and spatial location of 
small fires, and (2) we computed dNBR changes separately for 
FC in and FC out to estimate BA s f. With this method, we avoided 
having to design a detection algorithm that separated small 
fires from all other sources of surface reflectance variation. 
In effect, we were able to draw upon the skill of the MCD64A1 
burned area product in classifying large fires to simplify our 
analytical approach. Other approaches may be possible, 


including for example empirically relating FC out to high spa- 
tial resolution burned area observations, modeling burned 
area-fire size relationships using a combination of higher 
resolution Landsat imagery and coarser resolution burned area 
and land cover products, or by using fire radiative power. 
As discussed below, however, further reductions in uncer- 
tainty using any of these methods will be limited without 
new approaches for validation that specifically target regions 
where small fires are common. 

[ 59 ] Our study provides evidence that fires below the 
detection limit for MODIS MCD64A1 burned area can have 
important consequences for burned area and biomass burn- 
ing emissions. In North America, Portugal, and parts of 
Mali, by accounting for these fires, we improved the agree- 
ment between MODIS burned area and several aspects of 
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Year Year 

Figure 10. Burned area time series during 2001-2010 for BA MCD64Al (red diamonds), BA s f ( orange circles), 
and the combined total (blue triangles). The unit of measure is Mha/yr. 
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Figure 11 . Carbon emissions from fires time series during 2001-2010 for GFED3 simulations driven by 
BA mcd 64 a\ (red diamonds), BA s f (orange circles), and the combined total (blue triangles). The units are 
Tg C/yr. The regions are (a) global and (b) Southern Hemisphere South America, Southeast Asia, and 
Equatorial Asia. 


these regional burned area data sets. Specifically, for boreal 
forest fires in Alaska and Canada, our algorithm reduced the 
low bias observed in past work for MODIS burned area 
products [Kasischke et al., 2011] and slightly improved 


correlations describing the spatial pattern of burned area 
each year. For the continental U.S., significant improve- 
ments were observed for agricultural and prescribed fire 
burned area comparisons. In Portugal, the spatial correlation 


Table 4. Fire Radiative Power and Active Fires Within and Outside of MODIS Burned Areas 


Region 

Mean FRP 
for FC in 
(MW) 

Mean FRP 
for FC ou , 
(MW) 

Ratio 

(FRP ou ,/FRP ,-„) 

FC in 

(x 10 4 /Year) 

PC out 

(x 10 4 /Y ear) 

F C ou J F C to tal 

Cumulative 

FRPin 

(x 10 4 MW/Year) 

Cumulative 

FRP out 

(x 10 4 MW/Year) 

Ratio 

(Cum. FRP out / 
Cum. FRP in ) 

BONA 

50 

33 

0.65 

3.5 

1.4 

0.29 

176 

47 

0.26 

TENA 

48 

18 

0.37 

2.0 

5.0 

0.71 

97 

86 

0.89 

CEAM 

17 

12 

0.73 

1.4 

7.4 

0.84 

23 

92 

3.97 

NHSA 

16 

11 

0.71 

1.8 

5.0 

0.74 

28 

56 

1.99 

SHSA 

24 

15 

0.65 

15.7 

36.1 

0.70 

371 

557 

1.50 

EURO 

28 

13 

0.45 

0.7 

2.2 

0.75 

20 

28 

1.38 

MIDE 

20 

12 

0.62 

0.5 

1.7 

0.75 

11 

20 

1.87 

NHAF 

18 

11 

0.64 

42.5 

34.3 

0.45 

750 

386 

0.51 

SHAF 

20 

12 

0.62 

38.8 

35.5 

0.48 

755 

431 

0.57 

BOAS 

29 

19 

0.66 

8.3 

9.9 

0.54 

242 

189 

0.78 

CEAS 

28 

14 

0.52 

7.6 

15.6 

0.67 

212 

226 

1.07 

SEAS 

14 

12 

0.90 

3.3 

9.8 

0.75 

46 

120 

2.63 

EQAS 

16 

13 

0.81 

0.8 

3.6 

0.82 

12 

46 

3.73 

AUST 

29 

18 

0.61 

21.1 

5.9 

0.22 

615 

105 

0.17 

Global 

23 

14 

0.61 

148.1 

173.4 

0.54 

3357 

2388 

0.71 
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Table 5. Fire Carbon Emissions by Region Averaged During 
2001-2010 


Region 

C Emissions 
From BA MC D 64 Ai a 
(Tg C/yr) 

C Emissions 
From BA s f 
(Tg C/yr) 

C Emissions 
From BA tota i 
(Tg C/yr) 

Increase From 
Small Fire 
Emissions (%) 

BONA 

58 

5 

63 

9 

TENA 

10 

8 

18 

81 

CEAM 

24 

24 

48 

101 

NHSA 

28 

20 

49 

72 

SHSA 

267 

125 

391 

47 

EURO 

3 

3 

7 

98 

MIDE 

1 

1 

3 

90 

NHAF 

461 

124 

585 

27 

SHAF 

555 

154 

709 

28 

BOAS 

107 

38 

145 

35 

CEAS 

32 

35 

67 

107 

SEAS 

140 

66 

205 

47 

EQAS 

74 

41 

114 

55 

AUST 

122 

10 

132 

9 

Global 

1884 

654 

2538 

35 


a The carbon emissions for the baseline run (BA M cd64ai) were about 2% 
lower than the original GFED3 estimates from van der Werf et al. [20 1 0] at 
a global scale, although regional differences were larger and followed 
variations in burned area. 

between MODIS burned area and the national estimates 
increased in 7 of 8 years. In Mali, our small fire algorithm 
improved the agreement with Landsat-derived observations, 
with the largest relative adjustments occurring at the begin- 
ning and end of the fire season. More validation is needed, 
and observed differences between the spatial distribution of 
active fires and burned area may allow for the identification 
of regions for more intensive study. 

[ 60 ] The increases in carbon emissions associated with 
small fires may have implications for several carbon cycle, 
atmospheric chemistry, and human health issues. Higher 
levels of emissions are likely to reduce low biases in carbon 
monoxide and aerosols observed in atmospheric models 
driven by earlier versions of satellite-based fire emissions 
inventories [e.g., Tosca et al., 2012]. Also, fire contributions 
to variability in atmospheric CCD and CH 4 [e.g., Prentice 
et al., 2011] will likely increase when small fires are accu- 
rately classified and represented in biogeochemical models. 
Evidence for this comes from our GFED3 simulations that 
showed that the standard deviation of annual mean emissions 
increased in all regions except Northern Hemisphere Africa 
during 2001-2010 when we included BA s f in our simulations. 
In the tropics, the response of fire emissions to El Nino/ 
Southern Oscillation, the Atlantic Multidecadal Oscillation, 
and other climate modes [e.g., Chen et al., 201 1] is likely to 
be higher when small fires are integrated into emissions 
inventories because of the large relative increases in burned 
area in drought-sensitive forests and peatlands across South 
America and Southeast Asia (Table 3 and Figure 8). 

[61] Across different continents, carbon emissions increa- 
ses from small fires were higher in more densely populated 
areas. As a result, episodic and chronic exposure to fire- 
emitted particulate matter and ozone, and thus fire-induced 
human mortality [ Johnston et al., 2012], is likely to be con- 
siderably higher than for estimates derived from GFED3. 
An important next step in this context is to more systemati- 
cally quantify small fire impacts on atmospheric composition 
and exposure using a chemical transport model. Another 
critical issue is to assess overlap with inventories that include 


open field agricultural waste burning [e.g., Yevich and 
Logan, 2003]. 

4.2. Applicability to Other Burned 
Area Mapping Efforts 

[ 62 ] Do our results have applicability for moderate reso- 
lution burned area products other than the specific one used 
in our study? We believe that the degree of spatial agreement 
between active fires and burned area is a useful metric for 
validation of global burned area products, in addition to the 
more traditional approach of comparing moderate resolution 
burned area products with Landsat imagery and geographic 
information system bum perimeters. A key advantage of this 
metric is that global patterns can be assessed, and information 
about performance can be gained in regions where validation 
data are unavailable. Under ideal conditions, active fires 
should be within or near the perimeter of burned areas, taking 
into account the different spatial resolutions of the two pro- 
ducts. Divergence between the two products is expected 
when 1) the active fire product has false positive detections, 
2) thermal emissions from non-landscape fire sources have not 
been adequately screened, or 3) when fire size or fire severity 
is below the detection limit for the burned area algorithm. We 
believe that this third condition is the primary driver of the 
global patterns documented here in Figures 1 and 2, because 
false detections from MODIS active fires have been quantified 
and are mostly confined to specific landscape transitions 
[Schroeder et al., 2008] and because we attempted to carefully 
screen our database for non-fire hot spots (section 2.1). Fur- 
thermore, Schroeder et al. [2008] highlight a mechanism by 
which differential surface heating between vegetated and 
cleared areas can generate false detections, a problem that was 
more pronounced with Terra data than Aqua. In this study, the 
ratio of FC out to FC in was greater for Aqua than Terra. Other 
types of errors, such as omission errors with the active fire 
product or commission errors with burned area require diag- 
nosis using other approaches. Even with these limitations, we 
anticipate that future refinements to both burned area and 
active fire retrievals will lead to an improved consistency 
between these two types of data products, and thus the types of 
comparisons shown here may serve as a useful benchmark for 
product evaluation. 

[63] The spatial and temporal distribution of the fraction of 
active fires outside of burned area (Figures 1 and 2) provides 
some guidance for the design of future validation efforts. 
Important areas for better understanding small fire burned 
area include Guatemala, Mexico, India, Thailand, Indonesia, 
Ukraine, China, and the Southeastern U.S. These areas have 
not been targeted as extensively for field programs, with past 
efforts often focusing instead on savanna, boreal forest, and 
Mediterranean ecosystems where fires tend to be larger and 
to have more persistent effects on surface reflectance [ Justice 
et al., 1996; Loboda et al., 2007; Roy and Boschetti, 2009; 
Kasischke et al., 201 1], Successful field campaigns for small 
fire estimation will likely require a different investment 
strategy for project resources as compared with earlier 
efforts. For example, while Landsat imagery is likely to 
remain useful, a 16-day repeat cycle (at best) will not capture 
many bums in areas where the land is being intensively 
managed for crop or timber products. In this context, inter- 
national collaboration will be essential, with coordinated 
targeting of the same region with multiple commercial and 
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government sensors from several nations over a fire season. 
Frequent sampling with high spatial resolution imagery is 
also needed to resolve potential seasonal variations in the 
efficacy of burned area and active fire products caused by 
changes in cloud cover, phenology, aerosol loading, and 
surface soil moisture [e.g., Laris, 2005]. Parallel field sam- 
pling also is crucial for this class of fires, and requires 
ground-based measurements of bum perimeters, fire severity, 
and other forms of land management that might trigger false 
detections [e.g., McCarty et al., 2009], 

4.3. Comparison With Other Regional 
and Global Fire Studies 

[64] Several recent studies have combined active fire and 
burned area data sets to improve estimates of burned area and 
carbon emissions [Roy et al., 1999; Loboda et al., 2007; 
Boschetti and Roy, 2009; Giglio et al., 2009; Roberts et al, 
2011; Kaiser et al., 2012], Active fires have been used, for 
example, to refine thresholds for burned area detection [Roy 
et al., 1999; Loboda et al., 2007; Giglio et al., 2009], to esti- 
mate fuel consumption [Boschetti and Roy, 2009], and to more 
effectively estimate regional-scale biomass burning emissions 
[Roberts et al., 2011; Kaiser et al, 2012]. An important next 
step with respect to inter-comparing these different approaches 
is to reconcile our findings with the fusion approach proposed 
by Roberts et al. [2011]. In that study, the authors found that 
emissions for Africa increased by only 6% when SEVERI- 
derived fire radiative energy from fires detected outside of 
MODIS MCD45A1 burned area was included in their analysis 
system (they denoted as pathway 1). Here we used a different 
source for FRP and burned area, and estimate that 34-36% of 
MODIS-derived FRP was outside of burned areas for Africa 
(Table 4). A comprehensive analysis of the spatial distribution 
of MCD45A1 and MCD64A1 burned areas and their relation 
to MODIS and SEVER! active fire detections will be needed 
to understand these differences, building on work by Freeborn 
et al. [2009]. 

[65] In the context of other ongoing burned area and FRE 
studies, we believe our estimates of small fire impacts on 
global burned area and biomass burning emissions are likely 
conservative for four reasons. First, our algorithm reduced 
but did not eliminate low biases observed for MODIS burned 
area in several regions of North America and Africa. Second, 
some small fires are below the detection limit for the MODIS 
active fire product [Schroeder et al., 2008] and in our opin- 
ion, these fires are more likely to occur outside of existing 
burned area products derived from moderate resolution 
imagery. Third, our best estimate of global burned area and 
biomass burning emissions relied on observations from 
Terra. Use of early afternoon observations from Aqua con- 
siderably increased estimates of small fire burned area in our 
sensitivity study, primarily because proportionately more 
fires were detected outside of MCD64A1 burned areas. This 
is consistent with the observation that many agricultural and 
savanna fires tend to have a more pronounced diurnal cycle 
than extratropical forest fires, and one that peaks in early 
afternoon [Vermote et al., 2009; Mu et al, 201 1]. As a result, 
these fires may be more prevalent during the 1:30 pm Aqua 
overpass than during the 10:30 am Terra overpass. Fourth, as 
noted in the methods, we attributed all of the observed dNBR 
changes in equation (4) to burned area. If bum severity is 
lower (per nr of burned area) for small fires outside of 


MCD64A1 burned areas, then the small fire burned area 
fractions would have to be even higher to induce the same 
dNBR change. We note, however, that this latter scenario 
probably would not translate to higher carbon emissions 
because lower bum severities per unit of burned area would 
be expected to have lower levels of fuel consumption. 

4.4. Uncertainties 

[66] Our study reduced uncertainties in global biomass 
burning emissions estimates by improving the representation 
of a key class of fire sizes and by lowering burned area 
biases in North America and Africa. Nevertheless, uncer- 
tainties related to regional and global burned area estimates 
remain considerable and are related to multiple factors. For 
example, understory fires in tropical, temperate, and boreal 
forest biomes probably were not well represented by our 
approach and further work is needed to quantify these fires 
and integrate this information with burned area products for 
other fire types [Nepstad et al., 1999; Morton et al., 201 1]. 
Also, in many areas where biomass burning emissions are 
important, the seasonal dynamics of fuel properties and 
cloud cover and their subsequent impact on active fire 
detection efficiencies and the functioning of burned area 
algorithms remain poorly understood. 

[67] With the approach we developed here, important 
uncertainties were introduced by the spatial averaging of the 
500 m dNBR statistics to the ~l-km resolution of active fire 
detections, by our assumption of a linear relationship 
between dNBR changes and sub-pixel burned area, and also 
by the re-projection and resampling of the active fire swath 
observations on the MODIS sinusoidal grid that is an integral 
part of the MOD14A1 product. More specifically, the daily 
composition thermal anomaly/fire product (MOD14A1) may 
disproportionately increase FC out from resampling, but this 
bias also would likely decrease 7. Some of these errors could 
be avoided in a future reprocessing by using the swath-level 
fire location product that provides a more precise set of lati- 
tude and longitude coordinates for each active fire pixel 
(MCD14ML), and by filtering this product by a maximum 
threshold on the across-track scan angle. This latter constraint 
may reduce spatial smoothing effects that likely occur at high 
scan angles and also may provide a more representative set of 
active fire pixels for quantifying dNBR changes inside and 
outside of burned areas. By combining analysis of Landsat 
and MODIS in future work, it also may be possible to assess 
and modify the dNBR-fire size relationship, which is 
assumed to be constant here in our definition of 7. As a result 
of the comparisons presented here, additional modifications 
to the MCD64A1 burned area algorithm are currently being 
implemented to better represent fires in agriculture and areas 
of the tropics with dense cloud cover. 

[68] For our emissions estimates, important uncertainties 
exist with respect to the separation of the dNBR changes 
measured for FC out into burned area and fire severity 
(combustion completeness) components. To be conservative 
with respect to our estimates of small fire impacts on burned 
area, here we assumed the observed dNBR changes were 
related solely to burned area and as a consequence, we 
applied the same combustion completeness parameteriza- 
tions as in GFED3. While we do not believe this assumption 
led to a high bias in our carbon emissions estimates, addi- 
tional field observations of burned area perimeters, fuel 
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loads, and fuel consumption are needed to improve the 
parameterization of these processes in global models. 

5. Conclusions 

[69] We developed a unified approach for using informa- 
tion from active fire and 500 m burned area to estimate the 
contribution of small fires to global burned area and carbon 
emissions. We found good agreement between the spatial 
distribution of active fires and burned area in savannas and 
boreal forests. In other biomes, including croplands, wooded 
savannas, and tropical forests, we observed many active fires 
that were outside the perimeter of burned areas. We esti- 
mated the area burned for this set of small fires using a novel 
approach that required tracking the surface reflectance 
changes of individual active fire detections inside and out- 
side of burned areas, and in nearby control areas. At a global 
scale, accounting for small fires increased burned area and 
global carbon emissions by approximately 35%. Sensitivity 
analyses of key model parameters led to a range of estimates 
of global burned area increases from small fires, from 24% 
to 54%. Although preliminary, this analysis indicates that 
small fires are an important class of global biomass burning 
emissions and one that needs to be more specifically targeted 
in future research and validation efforts. 

[ 70 ] Multiple lines of evidence indicated that our estimates 
of small fire contributions to burned area and emissions were 
conservative. Uncertainties remain considerable for this class 
of fires, and our analysis provided some possible guidance 
with respect to the design of future validation programs. 
Specifically, to build reliable regional data sets for small 
fires, increased emphasis needs to be placed on the devel- 
opment of high temporal and spatial resolution satellite 
imagery time series (including acquisition and analysis of 
dozens, if not hundreds of images during a fire season in a 
specific study area) and simultaneous collection of fire peri- 
meters by mobile ground teams. Key directions for future 
research including quantifying the contribution of small fires 
to atmospheric trace gas variability, aerosol radiative forcing, 
and air quality. 
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